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ABSTRACT 

Correlations in the failures of subsystems or components may arise 
when they share a common, random environment, we snow' that a natural 
too! i or calculating the reliability of such systems is the characteristic 
’'jocf'ona of the random hazard rate. h(tj. Some general results for tne 
reliability of series and parallel systems in terms of the characteristic 
functional of the hazard rate are derived and applied to a number of 
models of random environments. The applications include r andom hazard 
functions arising from I) non-fatal shocks of random amplitudes, 2 ) a 
Markovi am Gamma-marginal stochastic process, 3) system $t r ess related 
tc un-repaued damage from incoming rounds, ana 4 t impulsive and 
acrumu'atea neat stresses from a laser battle. 



int.-cac. t icn 

Tonsider a system whose components are subject to a random 
environment 1 , possibly including random shocks, in general, the i- h 

component hazard rate, may be wmtt.en as the sum A *./<. part 
‘ r hp ‘ :•'$*. teeter^ mist ic) part. X f (t ). accounts 'or wearou* and those 
r 3nd ;m rv’ lures wmcn occur independently for each component. The 
second . stocnastic) part, qat,. arises mom those environmental 



condit'ons which are shared by two or more components in f he system 



we deveioD a general formalism for calculating the reliability of a 
system in a ranaom environment in terms of the characteristic 
functional of the hazard rate. Various models of random environments 
are then proposed and the corresponding expressions for the 
characteristic functionals are calculated. The reliability of systems in 
,uCh environments are then found using these characteristic functionals, 
me equivalent, deterministic hazard functions for single components in 
some of these models are found not to be constant in time 2 . 

In Section !i the general formalism is introduced and illustrated for a 
single component suD ject to an environment which includes non-fatal 
shocks with random amplitudes In Sections 111 and IV the reliability of 
series and parallel systems is presented in general terms and applied v o 
the ron-fatal snocK model. Section V presents a model for the hazard 
rate as a Markovian, Gamma-marginal stochastic process. Finally, in 
ItrCtion VI we present three models in which the failure rate is 
vporticrai to system stress. Each of these models can be applied to 
. ti-con''ponent systems using the methods developed in the earlier 
Sections. 

... One Component 

~ne simplest case is that of one component with hazard rate h(0, a 
non-negaiivs stocnastic process. Tne reliability of the component, 
subject to a particular realization of the environment (hence, h(t) ) is 
jef med as. 
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where T is the random variable (r.v.) equal to the time me component 
fails and the subscript denotes a particular realization of me stochastic 
process, h. Using the usual definition of the hazard rate 3 , this is. 

L 

R h = exp[ -J'h(t')dt' I. f?; 

o 

Hence, averaging over the random environments we have, 

t 

R(t) = E h [ R h ! - E{ exp[ -Jh(t')dt' 1 } (3) 

o 

or, 

R(t) = c h ln(oi | . ( J ) 



C n . t lT|i )! = E! expl i/T|(f) niriat ] ;, (5) 

0 

defines tne characteristic functional 4,5 of the hazard rate. In other 
words, this characteristic functional, evaluated for a particular value of 
the test Rjnct'on, r[(t), yields directly the reliability of a component. 

We will see below how this generalizes to multi-component systems in a 
random environment. 



We can write the reliability of a component or system in terms of an 
effective, deterministic hazard rate, h(t), using the usual definition of 
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hazard rate, 



hit) = - d/dt in { R(t) }, 

out we must remember not to use this "hazard rate" to calculate 
reliability of larger aggregations, unless such aggregations are composed 
of components or subsystems not snaring a common, random environment. 

Deterministic hazard rate Plus random shocks 

It' the environment gives rise to non-fatai random shocks with rate #(t), 
then the hazard rate may be represented by a deterministic function. 

Y.t;, rius delta function spikes, possibly with random amplitudes. For 
convenience we will sometimes replace A(t) or 'p(t) with constants. 
However, in all the reliability formulas derived here. At can always be 
rer deed with J X( b idf for non-constant background hazard rate, and 2ft 
can ne ''eplaced with j'dil )df for a non-homogeneous shock process. 

f.nat A l f i may be r eal!y an effective hazard rate, obtained from the 
preceding equation, when one portion of the random hazard rate is unique 
to eac n - q vidua I component, i.e. not shared among separate components. 

he razard rate for non-fatal shocks occurring at times t is, 
n 

nit) - a * “ exj S( t ’ - t (6) 
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where the pdf of the tj is 2f(t)/n , n - {2f(t‘)dt\ and n is a 

o 



Poisson r.v. with mean n 



The characteristic functional of this process, derived in Append’- A, s, 



Ch it [T((-)] = exp{ rx /T|(t')at' + J2s(t’j [C^C Tj(f') ) - Mat !, (7i 

O 0 



where C N is the ordinary character istic function of the shock 

amplitudes, . In particular, ;f the shock amplitudes have 3n 
exponential distribution (parameter 8) then, 

C ^( <*)) ~ 8/(8 - i ■!)], 
and 



t t 

C h *[T|( • :j - exp{ iX fri(t )dt + i/d(t') T|(t )/f8 - irjft )] dtx. 

0 0 

ui e r ce. setting T|(t‘) = i, trie reliability is. 



R(t) = exp{ -Xt - n /(8 + I] }, (10) 

or, assuming a stationary environment (25 - const.; and using & * 1 /o. 



R(t) - exp( -Xt - 2ft ca /[lx + i j ) 
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The effective hazard rate. h(t), is constant for this model. Note that as 
c* -> oo the rate of failure becomes X plus the rate, o, of occurrence of 
natal) shocks (in fact. X could be a background rate of infinitely ntgn 
shocks). 



II!. Two different components m f he same environment 
if two or more components are in different, independent environments 
the system reliability follows from the usual formulas 3 expressed in 
terms of the individual reliabilities. We consider here two different 
components seeing the same environment. The reliability of two 
components in series subject to a given realization of the environmental 
condemns 's. 

p W r -/ r P h^ 1 = R h,(0 ' R h 2 ^- 02) 

If the components are in a parallel, redundant system the reliability is, 



O = ’ - o - R h (t)! • [I - RkU)] 



03 ) 



exp; -j n , v t ;dt ] + exp[ -Jh 2 (t')dt ] - exp{ -J[h,(t') + h 2 ( t )] dt ’ ), 



w-'ere ail integrals are 0 -> t. Hence, averaging over realizations of the 
environment we have, using the definition of the characteristic 

: u not icna . , 



nv w 



(14) 






r, i 



i n (■ j- 1 



01 



•)! 






C -V 



■ !t|(-)! 



T)( )=l 



clearly, the last term also represents the reliability of a series circuit 
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with the same two components. In the following ;t will pe understood 
that all characteristic functionals are for the processes over the 
nterval [0, t]. 

As an example, consider two components in series, seeing the same 
shocks, Put experiencing different amplitudes, c*, k$; where c< and 3 are 
u.d. exponential r.v's. Then h, + h 2 is given by 

n 

h(t') = A| + X 2 + I (ry, * kj ) SCt* - l,). (15) 

1=1 

In this case C in Eq. (7) is replaced Py the characteristic function of 
the sum of two exponentials, 

C^vOl)) - o/[b " 10)] ' b / [ 0 " Kok), ( 16) 

and a is replaced by X, + X 2 , yielding, 

t 

C n , +h jri(-)] = exp{ i(X,+X 2 )jT|(t')dt‘ 1 • ( 1 7 1 

o 

t 

exp{ J 7>( t')dt { 8/(8 - rmt')] • 8/[8 - ’ rs ( t )kj - i] I. 
o 

*he re 1 1 at- 1 1 ' tu of the series system is then f 'ounc Py substituting 
7] ( t ) - :, yielding, 

R series (t) = ex P ( - ( VM f ' ^ 1 ’ & 2 /l(S + + ^ 1 )• (,?) 

The parallel system reliability is given Py, 
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G*> 



VO -- VO * R 2 (0 - Vries 



( 0 . 0 9 ) 

where R !>: are the same as Eq. (10) with S replaced with S and 8/k, 
respectively. 

if cx - 3 in the above, i.e. the amplitude of shocks seen by the two 
components are proportional, then we have, 

n 

h,(f) ♦ h 2 it‘) = X, + X 2 ♦ (Ok) 2 8(0 - t ( ). (20 J 

JO 

Tms is the same as for a single component, except X 0 _> X| + X 2 and 
-> s/( i *k). hence. 

VVV.)] r exp{ i(X/X 2 )jT|(t , )dt' } • (21) 

p 
j 

t 

L 

exp{ [2(t')dt'[ S/[8 - IT|(t‘) a: 0)1 - i 1 }, 

0 

and the "e ability of a series circuit is, 

V r!0S (O r exp{ - (X, + X 2 )t - 2ft [(k+!)/(8*k+1)] }. (22) 

The reliability of a parallel circuit is again given by Eq. (19). 

IV. m- Identical components 

For m components in series, subject to a given realization of tne 
environment, the reliability is, 
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( 25 ) 



z V° V° V 0 '" 



= exps -J(h, + • • • + h m )df }. 
o 

Hence, averaging over the environment, the rel.aDiiiLy is. 



R(t l 



-h, + lv 









( 24 ) 



a) !f all components respond to a given shock with independent 
amplitudes we obtain, 

R a (t) = exp{ -mX 0 t ♦ at [S m /(8 + !) rn - 1 3 }, (23) 



D; d the hazard rate is exactly the same for a 1 1 components, 

Rh/t) - expi -mX 0 t + 2t to/(o^rn ) - I] }. ’26) 



The effective hazard rate, h(t>. is again constant. Note that 
• I* 1/5) m = 1 + m/8 * • • •. 

nence, 

(i+ l/5) m > I + m/8, 

y % 

:l m /(S*l) m < S/(S + rn), 

which implies that R 3 < R t In other words, as expected, the reliability 

of a series system is higher in the more highly correlated environment. 
Also, using Eg. (19), it is clear that the reliability of a parallel system 
w > ’ i op (ewer in a more highly correlated environment. 
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m-identical components in parallel 

The reliaDiiity of m- identical components in a parallel redundant system 
subject to a given realization of the environment is, 

m t 

R h (t) = 1 - tt [ l - exD{-Jh j (t , )df} ]. (27) 

1= i o 

This can De expanded using binomial coefficients, , as, 
m 

R n U) = I - S(-)l C rri j ! R hi (t)”-R hJ U). (28) 

|=0 

Hence, averaging ove r the environment, we have, 
m 

wt; : - Z C m e~ iX o l exp( 2ft [SJ/(S*!)J - !] }. (29) 

1-0 

f or independent response to the shocks, and, 
m 

R(t) - I - Z (-) J C m , e'JV- exD( 3t [S/(S- j) - I] ), (30) 

J “0 

h there s exactly the same hazard rate for all components (each 
component sees exactly the same amplitude shocks). 

In general, tor m identical components in parallel, all with the same 
hazard rate, h(t) = X 0 (t) + g(t) (where g(t) is random), 
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m 



R(t) = I -Zi'-)lc m ) .e:<p!-ijx 0 (ndn C g [T|(')]| q( . )=r (31) 

1=0 

From the foregoing it is clear how to generalize the formulas for 
reliability in a random environment to more general configurations 
containing components both in senes and in parallel 

1 ) using the usual rules for probabilities, write the reliability of the 
system in terms of the individual component hazard processes, 
conditioned on a given realization of the environment, hence of the 
(h ] (t)} (cf. Eg. 27), 

2) average over the environment, hence over the [h^(t)j, 

3) express the resuit in terms of tne appropriate characteristic 
functionals (cf. Eq. 31), and finally, 

*4) obtain the effective component or system hazard rate, h(t). 

V. Expcnential/Gamma Hazard rate 

Consider an environment giving rise to a riarfov hazard rate process w 
a Gamma marginal ,,6 ’ r distribution which is common to all components 
.n tne system of interest. We show in Appendix B that the 
characteristic functional for such a process is, 

t. t t ' 

C_[qMl = explofdt' G(t"![3/[8 - 1 'dt."'g(t‘ )exp( - j m t ‘d c •] - I !) 

f 

I I 



o 



I. 



The correlation function for this process in the stationary case 
Co const.) is p r e"^ r . When u- 1 the process is exponential. 

T he reliability of a single component with hazard rate h(t) = X(t) + g(t) 
is then. 

R(f) = exp{ -JX(t')dt'} • Cg[T)(-)]| n(0 . (33) 

carrying out the resulting integrations yields. 

R(t) = (3d) 

expt - fX(t')df -0/>t/(2S + 1)} • {<58/(28*! - e"^ l )} u ^ 8+2)/(25S+ 1 ^ . 



T he effective hazard rate, h(t), is clearly not constant for this model. 
This can oe generalized to multi-component systems with the methods 
trom the previous sections. For example, the reliability of m identical 
proponents in series is obtained by replacing X -> mX and using 
m 1 - m-i in the characteristic functional. From Eg. (32) it is easily 
s^r-n that this leads to replacing 8 with 8/m in Eg. (34). 

c a 'yr 9 Rate proportional to System Stress 
c. this section we will consider three related models of hazard rate 
processes. In the first model we assume the rate of failure is 
proportional to the number of customers using the system, for example 
rn- rate oi wear on a highway may have a component which is 
proportional to the number of automobiles, N(t), using the highway. In 
the second model the failure rate of a major system (for example on a 
battleship during combat) is increased proportional to the amount of 
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un-repaired damage (the number of hits not repaired is N(t), the amount 
of aamage/hit is some positive r.v.). Both of these models are related 
to a generalized M/G/oo queue. Finally, we model the reliability of 
electronic systems in a laser battle scenario. Failures may be caused by 
impulsive stress caused by a laser hit, or by accumulated heat from the 
laser hits. It is found that none of the models in this section lead to a 
constant effective hazard rate. 



highway mode l 

Consider a large multi-lane highway in wmch the traffic level is such 
that cars do not interfere with eacn other. The number of cars on the 
highway can be modeled as an n/G/°o queue. The rate of arrivals of cars 
to the highway is 2f(t) and the pdf of the time spent on the highway 
section of interest is b(r) (the "service" distribution). The hazard rate 
is modelled as h(ti = X(t) + 6 N(t). We show in Appendix C that the 
characteristic functional of the hazard rate can be written as, 

t 

C n [T}(-)] = exp( -j X(t')T](t')dt'} • (35) 

o 

t oo t ' + x 

exp{ Jtf(t')dt'[l - JdcbUjexphjBr^sjds] I I. 
o o t 1 

For example, when the transit time is exponentially distributee, 
parameter u, the reliability can be evaluated as, 



R( fj 






xp( -20t/(p + 8j - tf£/(m 



3) 2 



-( n + 

- e c 



B)t; 



too; 
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And when the transit time is equal to x with probability one, the 
reliability is, 

R(t) = e"^ exp{ -3t - 3/0 • [i - e'^ r ] - 3(t-t)e"^ r }, (37) 

when t. > r, and 

R(t) = e" At exp{ -3t - 3/0 • [1 - e~^J } , (33) 

when t < r. These formulas generalize as shown above for components 
: n series and/or parallel systems. 

Battle Damage Model 

Assume that incoming rounds hit a ship according to a non-homogeneous 
Poisson process, rate 3(0. during a battle. The amount of damage to the 
ship per round (or, really, the increased system stress resulting from 
each hit) is given by a random variable, 0, with some non-negative 
distribution. Let the pdf of the time to repair the damage from each hit 
oe b ( 'T'. If the failure rate of the total system contains a term which 
increases proportional to the amount of un-repaireo damage (possibly 
because o* increased demands made on the rest of the system), then the 
"'azard r ate is, 

n 

hit ’ - X(t) ♦ k 2 0, {©( t - tj) - 0(t - t j -r,)}. (39) 



The stochastic portion of this corresponds, for example, to the total 
we gnt. of customers in a non-homogeneous f1/G/°° queue, where each 
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customer has a random weight. 8,. in Appendix C we derive the 

characteristic functional for this process. Applying that result we have, 
in general. 

t t ' + z 

R(t) = e'^ 1 exp{ -J’dtW) [ 1 - Cfj(jT|(s)ds)b(r)dr 1 1 1 ^ . (40) 

o t 1 

where Cg is the ordinary characteristic function of ft. If 8 is a constant 

this reduces to the previous model. Again, for m components in series 
replace X with mX and evaluate the characteristic functional for 
Tj (• ) - m-i. For parallel conf igurations Eq. (C4) can be used in 
conjunction with Eq. (31 ). 

One example that C3n be worked out in closed form is when 8 is an 
exponential random variable, parameter e. and the system cannot be 
repaired during the battle, i.e. b(r) = Sf z - <*>). The reliability of a 
Single component is easily calculated using Eq. (40) to be. 

R(t) = e~ ( ** +x)t - { (e - t)/e } eZ , (4! ) 

and for m identical components in series we obtain, 

R(t) = e '^ +mX ' jt • { (e - mt)/£ }^ /m . (42; 

_as>'-r -Satt m Model 

We now mode! the failures of electronic systems in a laser battle as. 
i > immediate failures resulting from impulsive thermal or f inetic 
shocks, or 2) random failures brought about by the accumulated thermal 
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stress from repeated hits by laser beams. The background hazard rate is 
again represented by X(t). Assuming an exponential cooling law, the 
hazard rate resulting from laser hits at times {tj} is, 

n 

hit) = Aft) * 2 {cxSU-tj) - 3 e(t - t,) e _x(t_t -j ) }, (43) 

r> 

where the second term represents the contribution to AT (increase in 
the system temperature) resulting from the laser strike at time tj. 

Using the same methods as in the Appendices it can be shown that the 
characteristic functional for the stochastic part of this process is, 

t t 

C[t; ■. • )! = exp{{dt'2f(t') [ exp[Jds T|(s) r(s - t')] - ’ll, ( 4 4) 

o r 

where, 

r ; s - t ) = c S(s-f) + £ e(s - t') e"^ 5 " 1 \ (45) 

Hence, the reliability is, 

= e"^* exp{Jdt*tf(t’) [ exp[ -cx - ft/'x-O - e' K ^ _! - Mi - l| ). (46) 

For m identical components in series, seeing exactly the same 
environment, Eg. (44) is modified by multiplying X, cx and $ by rn. The 
non-constant effective hazard rate, h(t), can be read directly from 
Eh. ’46! e g. for m identical components in series it is, 
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h(t) = mX(t) + 2f(t) { exp[ -mc< - iti0/k-( 1 - e ^ )1 - 1 }. (46) 



Note that d is not multiplied by m. All components see exactly the same 
shocks from laser strikes. 

More generally, we could model different types of laser weapons, ranges, 
atmospheric propagation, etc. by taking and $ as (correlated) random 
variables. The above equations can also be generalized if the 
components are not identical but have different vulnerabilities, thermal 
conductivities, etc. 

Vll- Conclusion 

We have presented a general formalism for calculating the reliability of 
multi-component systems subject to a random environment. The 
correlation in failures of different components can be accounted for by 
calculating the system reliability in terms of appropriate combinations 
of the characteristic functional of the random hazard rate, evaluated for 
T;(-) - i * j { j = 1 , n; i z /(-l)}. We have shown how the method may be 
used in deferent circumstances by applying it to a number of different 
models f or the random environment. Our results include non-constant 
effective hazard rates for some of the random environment models. Trie 
intuitively reasonable result that the reliability of series (parallel) 
systems is greater (less) in a correlated environment than when the 
components see independent environments has been demonstrated for 
some of our models. Although we have only considered series or parallel 
systems, the application to more complicated systems containing 
components both m series and in parallel is straightforward. 
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Appendix A. Derivation of the Characteristic Functional for Impulses 
we now derive the characteristic functional for the hazard rate resulting 
trom random delta function impulses. Unless otherwise noted, we define 
the character istic functional for the process over the interval C -> t. 
First, f— random part of the hazard rate, g(t), is defined more explicitly 
as, 

n 

gtf") - 3 cxj S(t' - tj), (A 1 ) 

.1 = 1 

where n has a non-homogeneous Poisson distribution, p n , with mean, 

f 

n - f 2f(t)dt, ; A '2 'i 

o 

the t, have pdf given Dy tf(C )/ n (as is appropriate tor a 

J 
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non-homogeneous Poisson process), and <x, have an arbitrary distribution, 
f(o< j ). with ordinary characteristic function. C^. The characteristic 
functional for g(t') is, 

t 

Cgl-qOl = E{ exp[ if^ft) g(t) dt] }, (A3) 

o 

or, oo n 

- 7 p n { fT [ JJdtj 2(t j)/n" dc<j f(c< j) e 1 *] ] ). 

n=0 j=1 

where the integration over tj is 0 -> t the integration over o<j is 0 -> <*>. 

We nave made use of the Dirac delta function to perform the integration 
over t\ The product over j now reduces simply to the expression in the 
sguare brackets raised to the rr h power because of the independence of 
each term in the snot noise- 1 ike process. Performing the average over n 
yielding the standard result for the generating function of a Poisson 

distribution) and the integration over <*, (yielding the characteristic 

function for the amplitude distribution) and taking the limit T --> oo, 
r 9suits in, 



t t 

C h r [ri( )] - exp{ iX 0 jT](t')df * {3(f) [C c< ( r\(V) ) - l]dt'}, (A4) 

O 0 

where wp have included the (independent) characteristic functional for 
tne deterministic portion of the hazard rate. 
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Appendix B. The Characteristic Function al of the Gamma Process 6 
Using the results of Appendix A, Eq. (Ad), the characteristic functional 
of random delta function impulses (rate oA(t) ) with exponential 
amplitudes (parameter 8). is. 

00 

Cg[r» (-)] = exp { i J dt oA(t) T) (t) / [8 - i T) ( t ) ] }. (Bi) 

“00 

If this shot noise is fed into a first-order linear filter (decay parameter 
A(t) ). the resulting process, x, satisfies, 

x + Ax = g. (B2) 

Now, to find the character istic functional for the process. x(t), subject 
to the initial condition, x(0) = x 0 , use the solution of the stochastic 
differential equation, (letting A be constant to make it easier to follow 
the derivation) in the definition, 

CO 

C x [rj(-)! = E{ exp[ i/n(t) x(t) dt] }, (B3 ■ 

o 

Or, 

= exp(ix 0 fdt Tj(t)e~^] • expiijdt rifUe'^-jdt 'e^ 1 j g(f') ). 

wnere we nave used the properties of the Heaviside unit step function, 
O(t-t'), so that all the integrations are over 0 -> <*>, thus making it 
easier to interchange the order of integration to give, 
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C X [T)(-)] = (B4) 

exp[ix 0 Jdt T((t)e“" t ] • exp{ i/dt g(t j [Jdte'^ t_t ) e(t-t') Ti(t)l }. 

The second exponential term is now of exactly the same form as the 
characteristic functional for g(t), with the expression in square brackets 
replacing Ti(t’). Hence, making use of Eq. (Bl) we have, 

00 

C sJtj (■)] = exp[ix 0 Jdt T](t)e"^ t -] * (B5) 

o 

00 00 00 

exp { i j dt ' vo [Jdt ^ rj 1 1)] / [S - i [Jdte"^"^ ^(t)] ] }, 

o f t’ 

where we have used the properties of the step function. 

To determine the marginal distribution of this process use 

■git) - T| q o( t _ t q ) , which recovers the ordinary characteristic function 

for xU 0 ). A straightforward calculation yields, 

T, flo j[r 10 ] - explrr^Aoe - ^^] { [3 - it | 0 e _2st o] / [s - it] 0 ] } u . (B6) 

~his is the characteristic function of a r.v. that, with probability 

■-a = has the value XQe’^ 1 - , and with probability (l-o<) is the sum of 

x 0 s -<5t and a r.v. with exponential distribution of parameter 8. 

Now, taking t 0 --> Eg. (B6) yields the marginal distribution of the 
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steady state process, 



c x(oo)tTl 0 ] - {S/ [8 - iTi 0 ]}° . (B7) 

i.e. the characteristic function of a Gamma distribution, as promised. 

A ppendix C. The Characteristic Functional of an M/G/co Queue, and 
Extensions 

The h/G/'eo queue may be modeled as shot noise impulses of unit 
magnitude 8 ( "customers" arrive with rate n is again Poisson ), 

n 

g(t) -- I SCt - tp, (CD 

1=1 

which have been filtered through a linear system with a random response 
function given by, 

r(t - t,J = e(t - tj) - e(t - tj - rr ( ). (C2) 

i.e. a unit height pulse with random duration, r ( (the "andom "service" 
time). Hence, for 3 ^ = 1, the number in the system at time t is. 



N(t) = _ 3j (e(t - tp - 9(t - tj -z j)). 



(C5) 



We have included the factor, so that we can allow 'he system 
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response to also have a random height, corresponding to a random weight 
(or damage, system stress, pollution, etc.) of each customer. Using the 
definition of the characteristic functional and following the same steps 
as in the previous Appendices, we obtain. 

oo n °° 

Z Pr, '■ n 1 J jfdt j2(t j)/FT d3 j f( 3 j) dzr,b(rj)exp[i3 j/dsT|(s) ©(s.tj.tj)] ). 
n=0 j=l o 

where ©(s.tj.rp = e(t - tj) - 0(t - tj -Zj). Again using the 
independence of each term in the product, and the fact that they are all 
the same, we can perform the summation over p n (Poisson), and the 

average over 3 (obtaining the ordinary characteristic function, ). 
Finally the characteristic functional may be written, 

t t'-r 

Irnt = exp{ -Jdt'tf(t') [I - C^(Jr[(s)ds)b(r)dr ] }. (C4) 

o t‘ 
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